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ABSTRACT 


A Comparison of Experimental and Theoretical Results 
for Leakage, Pressure Gradient, and Rotordynamic Coefficients 
for Tapered Annular Gas Seals. (December 1986) 

David Alan Elrod, B.S., Louisiana State University 

Co-Chairs of Advisory Committee: Dr. Dara Childs 

Dr. Clayton Nelson 

A brief review of current annular seal theory and a discussion 
of the predicted effect on stiffness of tapering the seal stator are 
presented. An outline of Nelson’s analytical-computational method 
for determining rotordynamic coefficients for annular compressible- 
flow seals is included. Modifications to increase the maximum rotor 
speed of an existing air-seal test apparatus at Texas A&M University 
are described. Experimental results, including leakage, entrance-loss 
coefficients, pressure distributions, and normalized rotordynamic 
coefficients, are presented for four convergent-tapered, smooth- 
rotor, smooth-stator seals. A comparison of the test results shows 
that an inlet-to-exit clearance ratio of 1.5 to 2.0 provides the 
maximum direct stiffness, a clearance ratio of 2.5 provides the 
greatest stability, and a clearance ratio of 1 .0 provides the least 
stability. The experimental results are compared to theoretical 
results from Nelson's analysis with good agreement. Test results for 
cross-coupled stiffness show less sensitivity to fluid prerotation 
than predicted. 
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NOMENCLATURE 

A,B - Fourier coefficients for rotor motion 

C,c - direct and cross-coupled damping coefficients (FT/L) 

e - surface roughness (L) 

Sq - displacement of seal rotor from centered position (L) 

f - fanning friction factor 

- inlet seal clearance (L) 

Hq - exit seal clearance (L) 

H - average seal clearance (L) 

K,k - direct and cross-coupled stiffness coefficients (F/L) 
k - entrance-loss coefficient 

M,m - direct and cross-coupled added-mass coefficients (M) 

• 

m « fluid mass flow rate (M/T) 
no, mo ■ Hirs' coefficients 

p » fluid pressure (F/L*) 

R - seal radius (L) 

R2 ' “ 2PU2C/U - nominal axial Reynolds number 

Rg = 2pUC/u - nominal Reynolds number 

U « / “ mean fluid flow velocity (L/T) 

U » mean fluid circumferential velocity (L/T) 

U2 - mean fluid axial velocity (L/T) 

X,Y ■ radial seal displacements (L) 
y - ratio of specific heats for air 
^o " ®o ^ ^r “ eQUilibrium eccentricity ratio 
M •» fluid viscosity (FT/L*) 


r 


- fluid density (M/L*) 

■ fluid shear stress (F/L*) 

- fluid preswirl angle (rad) 

- shaft angular velocity (1/T) 

- shaft precessional velocity (1/T) 

- fluid viscosity (FT/L*) 

Subscripts 


rotor 

stator 

axial 


circumferential 
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CHAPTER I 

INTRODUCTION 

Recent trends in turbomachinery design have been toward 
increased speeds, tighter clearances, and higher power to weight 
ratios. As a result, problems with vibration and dynamic stability 
have become more critical. To industrial plants, the cost of 
equipment outage can be devastating. The ability to predict and/or 
modify the dynamic characteristics of turbomachines can reduce this 
cost. 

Most vibrations detected in rotating machinery are caused by 
orbital motion of the rotor. If the frequency of this motion 
coincides with the rotational speed of the rotor, the vibration is 
called "synchronous". "Subsynchronous" vibration refers to an 
orbital frequency below the rotational speed. 

the most frequent cause of synchronous vibration is unbalance in 
the rotor. The response to unbalance is also called forced vibration 
because the unbalance exerts a rotating force on the rotor as it 
spins. Synchronous vibration amplitude reaches a relative maximum 
when the rotor rotational speed coincides with a damped natural 
frequency of the rotor. Damping in the rotor system helps to reduce 
the peak synchronous vibration levels. If critical speeds must be 
traversed during start-up and shut-down of a turjjomachine, the 
designer aspires to introduce additional damping into the system. 
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Subsynchronous vibration problems are far less frequent in 
practice than synchronous. Self-excited motion is a type of 
subsynchronous vibration which can be particularly destructive. It 
is often the result of a tangential force acting on the rotor due to 
some fluid or friction mechanism. This motion occurs at a natural 
frequency of the rotor below the running speed. Motion which is 
self-excited appears suddenly at some threshold speed with large 
amplitudes which sustain or grow as the running speed increases. The 
consequences range from the prevention of running speed increases to 
catastrophic failure of the equipment. 

The forces developed in annular seals can play a role in both 
synchronous and self-excited vibration. Black et al. [1 - 3] have 

demonstrated the critical effects that forces developed by neck-ring 
and interstage seals have on the rotordynamic behavior of pumps. 
Also, stability difficulties with the Space Shuttle Main Engine 
(SSME) high-pressure fuel turbopump [^] have prompted further 
research into the forces developed by liquid seals. 

Experiences have shown that various gas seal configurations can. 
similarly influence the rotordynamic behavior of turbomachines. For 
example, in the high-pressure oxygen turbopump (HPOTP) of the SSME, 
initial vibration problems were remedied by changing the turbine 
interstage seal from a stepped-labyrinth configuration to a 
convergent taper seal with a honeycomb stator and a smooth rotor [5]. 
Compressible flow in annular seals has been analyzed by Fleming [6, 
7] and Nelson [8 - 10]. Both Fleming [6] and Nelson [9] suggest that 


the Inlet to exit clearance ratio In a tapered seal can be optimized 
to improve rotordynamic behavior by maximizing seal stiffness. 
Nelson et al. [11] have reported experimental leakage, stiffness, and 
damping data on two smooth rotor/smooth stator seals in which a 
tapered seal with an inlet to exit clearance ratio of 1 .5 was shown 
to be stiffen than a straight bore seal. No dependence of seal 
behavior on rotor speed was noted through the 8000 rpm maximum speed 
of these tests. Comparisons of the experimental data to the 
predictions of Nelson's analysis [9] showed Nelson underpredicts 
stiffness. Furthermore, the underprediction grew worse as the 
relative roughness, e/2HQ, of the seal increased. There have been no 
comparisons of experimental data to Nelson's most recent analysis 
[ 10 ]. 

The purpose of this report is threefold. Modifications on an 
existing air seal test apparatus at Texas A & M University are 
described, which make seal tests at rotor speeds up to 16000 rpm 
possible. A comparison of experimental data is provided from the 
tests of four smooth rotor /smooth stator seals with different inlet 
to exit clearance ratios. A comparison is provided between 
theoretically-predicted and experimentally-obtained data for the four 
seals. The leakage of the working fluid through the seal, the 
pressure gradient along the seal length, entrance pressure-loss data, 
and rotordynamic coefficients provide a basis for comparison. A 
short discussion on seal theory is included, as is a description of 
the rotordynamic identification scheme used. The work presented 


herein is intended to add to the database on seal forces and to 
determine the validity of a theoretical analysis for predicting those 


forces. 



CHAPTER II 


ANNULAR SEAL ANALYSIS REVIEW 


As related to rotordynamics, seal analysis has the objective of 
determining the reaction forces acting on the rotor arising from 
shaft motion within the seal. Due to similarities between plain 
Journal bearings and annular seals, seal analysis is generally based 
on governing equations which have previously been developed for 
bearings. 

Annular seals and plain bearings are geometrically similar, but 
seals typically have radial clearance-to-radius ratios on the order 
of 0.005, versus Cp/R ratios of 0.001 for bearings. Due to seal 
clearances and pressure differentials, fully-developed turbulent flow 
normally exists. Also, seals are nominally designed to operate in a 
centered position. Journal bearings, on the other hand, have 
operating eccentricities which vary with running speed and load. 
Therefore, most of the rotordynamic work for bearings has been done 
to determine dynamic coefficient versus eccentricity relationships. 

Two linearized seal models, expressed in terms of dynamic 
coefficients, have been suggested for the motion/reaction-force 
relationship. For small motions of the rotor about an arbitrary 
position in the seal, as shown in Fig.l, the relation can be written 
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Fig. 1 Sull •otlon of « aeal rotor about an aocantrlc poaltlon 
(m la the rotor apln apaed. Q ia the preceaalonal orbit 
frequency). 



Fig. 2 Snail ■otlcn of a acal rotor about a oantarad poaition 
(« la the rotor apln apeed, 0 ia the praoaaalonal orbit 
frequency ) . 



In this equation, {X, Y} define the motion of the seal rotor relative 
to its stator, {Fx» Fy} are the components of the reaction force 
acting on the rotor, and the dynamic coefficients (Kyy* Kyyi Cyx* 
Cyy, Mxx* Myy} and iKxyt Kyx» Cxy» Cyx* Mxy, Myx) represent the 
"direct" and "cross-coupled" stiffness, damping, and added-mass 
terms, respectively. These coefficients are functions of the 
equilibrium eccentricity ratio Cq = Sq / Cp, where the eccentricity 
ratio Eq equals the displacement (cq) of the rotor from the centered 
position divided by the nominal radial clearance (Cp). The term 
"cross-coupled" refers to the coupling effect exhibited by the 
off-diagonal terms; specifically, motion in one plane introduces 
reaction forces in an orthogonal one. These cross-coupled terms 
arise from the fluid's circumferential velocity component, and show a 
strong dependency on both the magnitude and direction (with respect 
to rotor rotation) of the fluid velocity. This circumferential 
velocity component may arise from the prerotation of the fluid as it 
enters the seal due to some rotating element upstream, or it may 
develop as the fluid passes through the seal, with rotor shear forces, 
"dragging" the viscous fluid around its periphery. The cross-coupled 
stiffness term usually produces a destabilizing force component, and, 
therefore, is of considerable interest. The cross-coupled damping 
and added-mass terms are generally much less influential than the 
cross-coupled stiffness term with respect to stability. For no fluid 
rotation, these cross-coupled terms are zero. 

The second linearized seal model applies for small motions of 


the rotor about a centered position in the seal, as shown in Fig. 2. 
This model can be expressed 



where the dynamic coefficient matrices are skew-symmetric. 

Theoretical work on annular seals has been done for both 
incompressible eind compressible fluids. Black et al. [12] have 
developed analytical "short-seal" solutions for incompressible seals, 
which account for circumferential fluid flow due to wall shear 
stresses but not pressure perturbations. The analysis employs a 
bulk-flow assumption and accounts for fluid prerotation as it enters 
the seal. Childs' [13] Incompressible seal analysis provides 
"finite-length" solutions, in which both shear and pressure-induced 
flow are included. Childs utilizes Hirs' [1^] turbulent bulk-flow 
model, and accounts for inlet swirl as well as perturbations in axial 
and circumferential Reynolds numbers due to clearance perturbations. 

Compressible flow in seals has been analyzed by Fleming [6, 7] 
and Nelson [9] . Fleming presents a short seal solution for the 
leakage, direct stiffness, and direct damping coefficients of 
straight and tapered, smooth, annular gas seals. However, he does 
not include the effects of the circumferential momentum equation, and. 
thus cannot obtain the cross-coupled terms. Nelson includes the 
effect of circumferential momentum in analyzing both smooth and 
surface-roughened annular seals in the straight and tapered 


configurations, and he solves for the cross-coupled coefficients. 
Nelson has developed two analyses, one of which [8, 9] is based on 
Hirs' [14] turbulent bulk-flow model. The other [10], which is used 
for comparison in this report, is a bulk-flow analysis which utilizes 
Moody's friction factor equation [15]. An outline of Nelson's 
analysis is included in the section that follows. 

One interesting result of Fleming's analysis which was later 
supported by Nelson's analysis concerned the behavior of convergent- 
tapered seals. A convergent-tapered seal is one which has an inlet 
clearance which is larger than the exit clearance. Both Fleming and 
Nelson predict that if the inlet clearance of a constant clearance 
annular seal is increased while the exit clearance is held constant, 
the direct stiffness of the seal will increase to a point and then 
decrease as the inlet clearance is increased further. Nelson et al. 
[11] have reported on the experimental verification of an increase of 
stiffness with taper (inlet to exit clearance ratio). An 
experimental investigation of the effect of taper on rotordynamic 
coefficients is included in this report. 


CHAPTER III 


NELSON’S ANALYSIS 

Nelson [8-10] has developed analyses which provide both static 
and dynamic results for annular gas seals. The static results 
Include fluid leakage through the seal, pressure gradient along the 
seal axis, and the fluid axial and circumferential velocities through 
the seal. Dynamic data provided by the analyses consist of the 
rotordynamic coefficients (direct and cross-coupled stiffness and 
damping terms) for small rotor motion about a centered position 
(equation (2)). Nelson assumes that the added-mass terms are 
negligible for gas seals, and, hence, equation (2) is written 



In one approach, Nelson [8,9] utilizes a modified Hirs' [1M] 
turbulent bulk-flow fluid model to develop • governing axial and 
circumferential momentum equations, and his model is completed by the 
continuity and energy equations. Hirs' model assumes that the rotor 
and stator friction factors (fp and fg, respectively) can be written 
as; 

f - no(Re)“o (^) 
where Rg is the Reynolds number relative to the surface upon which 
the shear stress is acting, and the constants no and mo are generally 
empirically determined from static pressure flow experiments. 


Comparisons of the predictions of this model with experimental 
results [16] show that the model underpredicts the experimentally 
measured direct stiffness. Furthermore, this underprediction worsens 
as the relative roughness (e/2H) of the seal increases. 

In another bulk-flow analysis. Nelson [10] uses Moody's pipe 
friction formula [15]. Previously used by von Pragenau [17] in an 
incompressible seal analysis, this formula approximates the Moody 
chart: 

6 1/3 

f - 0.001375 { 1 > [20000 e/2H + (10) /Rg] 1 (5) 

where e is the roughness ("bump" height), H is the seal clearance, 2H 
is the hydraulic diameter, and Rq is the Reynolds number relative to 
the surface upon which the shear stress is acting. This equation 
gives friction factors within 5% of those in Moody's diagram for 
Reynolds numbers between 4000 and 10’ and values of e/2H up to 0.01. 
For e/2H > 0.01, equation (5) significantly underestimates the 
friction factor. The Reynolds numbers at the stator, Rgg, and rotor, 
Rgp, are 

2 2 1/2 

Res - P (2H) (Uz + Ug ) / y 

2 2 1/2 

Rgr “ P (2H) [Uz + (Uq - Ru) ] / u 

where Uz and Ug are the axial and circumferential components of the 
mean fluid velocity, p is the fluid density, y. is the fluid 
viscosity, R is the seal radius, and u is the rotor speed in radians 
per second. As applied to an annular pressure seal, fg and fp 
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For this report, the rotor and stator surface roughnesses were 
determined using a profilometer. 

Assuming small motion of the rotor about a centered position. 
Nelson uses a perturbation analysis similar to that employed by 
Childs [13] to develop zeroth ^d first-order perturbation equations. 
The zeroth-order solution represents a zero-eccentricity flow 
condition, with rotor rotation but without precession. This solution 
is developed iteratively and yields the mass-leakage flow rate, and 
the axial distribution of pressure, axial velocity, density, and 
circumferential velocity. The iterative solution scheme uses initial 
guesses for the zeroth-order seal entrance Mach number and entrance 
pressure-loss coefficient. The entrance-loss relationship is defined 
by 


V(1-Y) 

Po(0) - {1 + C(T-1 )(k+1 )Mo2(0)] / 2} (8) 


where Po(0) is the seal entrance/reservoir pressure ratio and Mo(0) 
is the entrance Mach number. The entrance Mach number is iteratively 


adjusted, and the loss coefficient k is recalculated for a given seal 
according to the empirical relationship 

k - a ♦ b(Rz)C (9) 

In equation (9), Rz is the axial Reynolds number, and a, b, and c are 
constants calculated from test data and provided as input paraimeters 
for Nelson's ainalysis. The Iterative solution procedure for Mq( 0) and 
k continues until either: 

1 ) the Mach number at the exit reaches unity and the exit 
pressure is greater than the sump pressure (choked flow) , or 

2) the exit pressure equals the sump pressure and the exit Mach 
number is less than unity (unchoked flow). 

The pressure, density, and velocity distribution and their 
derivatives which are determined in, the zeroth-order solution, and 
the entrance-loss relationship determined from test data, are used in 
defining coefficients of the first-order perturbation equations. 
These equations define the pressure, density, and axial and 
circumferential velocity perturbations due to rotor motion. The four 
physical boundary conditions required for the solution of these 
equations depend on the perturbation conditions that are specified at 
the seal entrance and exit. The inlet circumferential velocity 
perturbation is zero. Expansion of the entrance pressure-loss 
relationship of equation (8) yields a second boundary condition. A 
similar expansion of an equation for the density change at the 
entrance provides another boundary condition. Finally, for choked 
flow, the first-order perturbation in the exit Mach number is zero. 


while for unchoked flow, the first-order perturbation in the exit 
pressure is zero. 

Application of these boundary conditions and numerical 
Integration of the ordinary differential equations provides the 
first-order solution. Integration of the first-order pressure 
solution along and around the seal periphery yields the direct and 
cross-coupled stiffness and damping coefficients, K, k, C, and c, 
respectively. 

The input parameters which can be varied in Nelson's analysis 
Include: 

1) reservoir pressure and temperature, 

2) sump pressure, 

3) seal geometry (i.e. radius, length, clearances), 

M) rotor rotational speed and precession rate, 

5) entrance circumferential velocity of the fluid, 

6) rotor and stator surface roughness, 

7) empirical entrance-loss relationship,. 

8) fluid viscosity, gas constant, and ratio of specific heats. 

An evaluation of Nelson's new model using Moody's equation can 

be accomplished using the Texas A & M University gas seal test 
apparatus described by Childs et al. [18]. The apparatus has 
previously been used in the evaluation of Nelson's Hirs-based model 
[11, 16], as well as Scharrer's Iwatsubo-based model for labyrinth 
seals [19]. This report includes an extensive comparison of the 
predictions of Nelson's new model to test results obtained using the 


TAMU apparatus and one example comparing the predictions of 
two analyses. 


Nelson's 
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CHAPTER IV 

TEST METHOD AND APPARATUS 

Test Method. The test method used In this study for determining the 
rotordynamic coefficients of seals has been used previously by lino 
and Kaneko [20] . An external hydraulic shaker is used to impart 
translatory harmonic motion to the rotating seal, and rotor motion 
relative to the stator and the reaction force components acting on 
the stator are measured. 

Fig. 3 shows the manner in which the rotor could be positioned 
and oscillated in order to identify the dynamic coefficients of the 
seal for small motion about Cq. If the added-mass terms are assumed 
negligible, equation (1 ) is rewritten 

+ . ( 10 ) 
*^Yx(eo^ Sy(':o)_ X Sx^E^o) ^YY^^o^J ^ 

First, harmonic horizontal motion of the rotor is assumed, where 

X = eg + A sin(Qt) + B cos(Dt) 

X - AQ cos(Qt) - Bn sin (nt) 

Y - Y - 0 . 

This yields small motion parallel to the static eccentricity vector, 
where Q is the shaking frequency. In a similar fashion, the X and 
Y-direction force components can be expressed 

^X “ ^Xo ^XS sin(nt) + ^XC cbs(Qt) 

Fy ■» ^Yo ^YS sin(nt) + ^YC cos(nt) . 



( 11 ) 



Substituting these expressions into equation (10) and equating 
coefficients of constant, sine, and cosine terms yields the following 
four equations for the dynamic coefficients 

■ fxS ” Kxx A - Cxx B Q 

■ f'xc “ Kxx ® Cxx A Q 

( 12 ) 

- Fys “ Kyx A - Cyx B n 

- Fyc “ Kyx B + Cyx A 0 . 

Solving this system of four equations in four unknowns defines the 
dynamic coefficients as 

*^Xx(^o) “ “ ^^XC B + ^XS A) / (A* +. B*) 

Sx(eo) - - ^^YS A + ^yc B) / (A* + B*) 

(13) 

^XX^^O^ “ ^XS B - ^XC A) / Q(A* + B*) 

^YX^^O^ ” ^^YS ® ~ ^YC A) / D(A* + B*). 

Therefore, by measuring the reaction forces due to known rotor 
motion, determining the Fourier coefficients (A, B, Fys* Fxc> ^YS* 
Fyc), and substituting into the above definitions, the indicated 
dynamic coefficients can be. identified. If the rotor is shaken about 
a centered position (sq » 0), then the process is complete. Since 
the linearized model has skew-symmetric stiffness and damping 
matrices, all of the coefficients are identified. If, however, the 
rotor is shaken about an eccentric position, as initially postulated, 
then it must be shaken vertically about that same point in order to 
complete the identification process. 

Assuming harmonic vertical motion of the rotor, as defined by 
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X - Cq, X - 0, 

Y - A 3ln(Qt) ♦ B cos(Qt), 

Y - AQ eos(Qt) - BQ sin(Qt), 

yields oscillatory motion that is perpendicular to the assumed static 
eccentricity vector. A similar process as before results in the 
coefficient definitions 

*^YY^^o) “ " (^YS * * ^YC ^ + B*) 

®S(Y^*=o) - - (^XS A ^XC ^ (A* 

(U) 

^YY^^O^ “ (^YS ® ~ ^YC ^ Q(A* + B*) 

^XY^^o) ” (^XS ® “ ^XC A) / Q(A* + B*) . 

All eight dynamic coefficients are thus determined by alternately 
shaking the rotor at one frequency Q in directions which are parallel 
and perpendicular to the static eccentricity vector. 

Test Apparatus. The test apparatus used for this study was the Texas 
A & M University gas seal apparatus described by Childs et al. [18] 
modified for rotor speeds up to 16,000 cpm. Fig. M shows the 
assembled rig. The test apparatus allows the following seal 
parameters to be controlled independently, providing insight into the 
influence these parameters have on seal behavior: 

1) static position and dynamic motion of the rotor, 

2) pressure ratio across the seal, 

3) prerotation of the incoming fluid (air), 

. • 

4) seal configuration, and 

5) rotor rotational speed. 

Previously, the maximum possible test speed was 8500 cpm. High 
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bearing temperatures and the reduction of interference in the rotor- 
shaft fitment with increasing speed had served to limit shaft speed. 
These problems have been addressed by some specific design 
modifications which are discussed below. 

In the past, the seal rotor was press-fitted and secured axially 
by a bolt circle to the rotor shaft. As the running speed is 
increased, however, the inertia-induced diametrical growth of the 
rotor exceeds the growth of the shaft. By increasing the interference 
in the stationary rotor-shaft fit, a greater allowance for this 
growth difference has been provided. Fig. 5 shows the present rotor- 
shaft design, a tapered rotor which is hydraulically expanded during 
installation. The rotor is inserted over the end of the tapered 
shaft and a large nut is used to pull the rotor onto the shaft. 
Fluid is pumped between the shaft and rotor, causing the rotor to 
expand. This separating force allows the rotor to be pulled onto the 
shaft until the desired interference fit is achieved. 

The problem of high bearing temperatures has been eliminated by 
replacing a roller-type thrust bearing and modifying the lubricant 
flow. A Torrington Hydraflex thrust bearing, consisting of eight 
one-inch rubber-faced pads which are water lubricated, is now in 
place at the rear of the rotor. In addition, the lubricant for the 
Torrington hollow-roller bearings which support the shaft has been 
changed to light turbine oil with a maximum temperature of 270®F. 
The hollow-roller-bearing caps have been modified to direct the oil 
flow to the regions of heat buildup. These modifications are shown 


HydrftttUcttlty Fllttil 
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Cross-sectional view of test section 
showing rotor-shaft assembly. 




The final modification to allow operation of the TAMU gas seal 
test apparatus at high speeds was the installation of Koppers 
circumferential seals for the hollow-roller and thrust bearing 
lubrication systems. At 16,000 rpm, the surface speeds of the shaft 
and rotor (170 and 350 ft/sec, respectively) exceed the limits of lip 
seals, which had been used on the TAMU apparatus. The Koppers seals 
in Fig. 5 were designed for gas applications. The sealing mechanism 
is a segmented carbon seal ring. 


CHAPTER V 


TEST RESULTS: INTRODUCTION 

The results reported here are from tests of four smooth- 
rotor/smooth-stator seals with different ratios of inlet-to-exit 
clearance (i.e. different tapers): 1.0, 1.5, 2.0, and 2.5. For the 

remainder of this report, the seals will be referred to as seal 1.0, 
seal 1.5, seal 2.0, and seal 2.5. Fig. 6 shows a tapered stator. 
Tables 1 - 3 contain the dimensions and pertinent data for each seal 

tested. 

The test program had the following objectives: 

1) Acquire leakage, stiffness, and damping coefficients as a 
function of rotor speed, pressure drop, and inlet circumferential 
velocity for four smooth-rotor /smooth-stator seals with different 
tapers but equal exit clearances. 

2) Compare the effect of varying the seal taper angle on the 
experimentally determined rotordynamic coefficients. 

3) Compare test results to the predictions of Nelson's analysis 
which uses Moody's equation for a friction factor. 

When shaking about the centered position, the test apparatus can 
be used to control the rotor speed, reservoir pressure (i.e. supply 
pressure), circumferential velocity of the inlet air, and the 
frequency and amplitude of translatory rotor motion.' Two shake 
frequencies, 58.8 and 7^.6 Hz, were used during testing with 
essentially the same results. The results plotted here were obtained 







Table 1. Teat atator apecificationa. 


Statora 

Seal: 


1.0 

1 .5 

Diameter: 

upatream 

15.22 

cm (5.991 in) 

15.26 cm (6.008 in) 

downatream 

15.22 

cm (5.991 in) 

15.22 cm (5.992 in) 

Material : 


braaa 

brass 

Surface 
roughneaa ; 

0.406 

i pm (16 pin) 

0.660 pm (26 pin) 

Seal: 


2.0 

2.5 

Diameter: 

upatream 

15.29 

cm (6.020 in) 

15.34 cm (6.040 in) 

downatream 

15.22 

cm (5.990 in) 

15.22 cm (5.992 in) 

Material: ' 

. 

brass 

brass' ~ 

Surface 

roughneaa: 

0.381 

pm (15 pin) 

0.559 pro (22 pin) 


Table 2. Teat rotor apecificationa. 
Seal; 1.0, 1.5, 2.0, 2.5 

Diameter : 

upatream 15.14 cm (5.959 in) 

downatream 15.14 cm (5.959 in) 


Material: 304 atainless steel 

Surface 

roughneaa: 0.279, pm (11 pin) 



Table 3. Test seal specifications. 


Seal 

Seal 

length 



Radial 

clearance 

Hj^ 

1 .0 

5.080 

cm 

(2.000 

in) 

Hi 

m 

0.406 

mm 

(0.016 in) 

1 .0 






Ho 

- 

0.406 

mm 

(0.016 in) 


1 .5 

5.080 

cm 

(2.000 

in) 

Hi 

. 

0.622 

mm 

(0.0245 in) 

1 .48 






Ho 

a 

0.419 

mm 

(0.0165 in) 


2.0 

5.080 

cm 

(2.000 

in) 

Hi 

m 

0.775 

mm 

(0.0305 in) 

1 .97 






Ho 


0.394 

mm 

(0.0155 in) 


2.5 

5.080 

cm 

(2.000 

in) 

Hi 

a 

1 .030 

mm 

(0.0405 in) 

2.45 






Ho 

a 

0.419 

mm 

(0.0165 in) 



Table 4. Definition of symbols used in figures. 


Pressures 

1 - 3.08 bar 

2 - 4.46 bar 

3 - 5.84 bar 

4 - 7.22 bar 

5 - 8.25 bar 


Rotor speeds * 

1 - 3000 cpm 

2 - 6000 cpm 

3 - 9500 cpm 

4 - 1 3000 cpm 

5 - 16000 cpm 


Inlet circumferential velocities 

1 - High velocity against rotation 

2 - Low velocity against rotation 

3 - Zero circumferential velocity 

4 - Low velocity with rotation 

5 - High velocity with rotation 


The pressure for each test is set at the flowmeter of Fig. 7. 


* For high inlet circumferential air velocity against rotation, seal 
1.5 was tested at the following rotor speeds: 3000, 6000, 8000, 
10000, 12000, 14000, and 16000 cpm. 



by shaking at 7^.6 Hz at an amplitude between 4 and 5 mils. The 
actual test points for each of the other three independent variables 
are shown in Table 4. 

Fig. 7-10 show the inlet circumferential velocity values (U.) 
for the configurations described in Table 4 for the seals reported on 
here. The equation for U. is 

Ug - m sin 4> / p Ay 

where m is the fluid mass flow rate, p is the fluid density. Ay is 
the exit area of the fluid turning vanes, and 4) is the fluid swirl 
angle at the turning vanes exit measured from the axial direction. 
Negative circumferential velocities represent velocities opposed to 
the direction of rotor rotation. Positive velocities are in the 
direction of rotor rotation. Note that curve 3 (representing zero 
inlet circumferential velocity) lies on the horizontal axis in each 
figure. The inlet circumferential velocity ratio, the ratio of inlet 
circumferential velocity to rotor surface velocity, ranged from about 
-6 to about 8. When reviewing the following figures. Table 4 and 
Fig. 7-10 should be consulted for the definitions of symbols used. 
Normalized Parameters. Before the tests described herein were 
performed, the TAMU gas seal test apparatus was modified as described 
in the TEST METHOD AND APPARATUS chapter to allow operation at 
running speeds up to 16,000 cpm. As expected, subsequent tests 
revealed a dependence of the rotor diameter on running speed due to 
inertia and thermal effects. The rotor growth data are shown in 
Table 5. Thus, as the rotor turns faster, the forces in the seal are 


CmCUMFERENTIAL VELOC-ITY CrRCUMFC^»CM^f AL VELOCIfr (m/aae) 



PRESSURE RATIO <bor/bor> 

Fig. 7 Inlet circumferential velocity vs. pressure 
ratio for straight bore seal (seal 1.0) 
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Table 5. Q'owth of rotor with rotational speed. 


Rotor speed 
(rpm) 


Diametrical growth 
(mm) (inches x 1000) 


3000 

6000 

9500 

13000 

16000 


0.01 

O.M 

0.02 

0.7 

0.03 

1 .2 

0.05 

1.9 

0.11 

M.4 


Table 6.. Normalized coefficients. 

Ho 

K * K , (nondim) 

DL(AP) 

Ho 

C » C (sec) 

DL(AP) 

“ stiffness (N/mm) C - damping (N sec/mm) 

» seal exit clearance (mm) D - seal diameter (m) 

- seal length (m) AP « pressure drop across 

seal (N/m*) 


affected not only by the increased surface speed of the rotor (drag) 
but also by a. change in clearance (friction factor). The 
rotordynamic coefficients presented in this study have been 
normalized. See Table 6 for the definitions of the normalized 
parameters. Note that the effect on the physical rotordynamic 
coefficients of increasing the running speed is generally greater 
than on the normalized coefficients. Fig. 11 shows the dependence on 
running speed of the direct stiffness of the 1 .0 taper seal before 
and after normalization. Note the smaller percentage increase of the 
dimensionless stiffness as the speed increases. 

Relative uncertainty. Before the test results are given, a statement 
about the experimental uncertainty is needed. The method used is that 
described by Holman [21] for estimating the uncertainty in a 
calculated result based on the uncertainties in primary measurements. 
The uncertainty wr in a result R which is a function of n primary 
measurements xi ,X 2 ,X 3 , . . . ,Xp with uncertainties wi ,W 2 ,W 3 , . . . .w^ is 



1/2 


(15) 


In this case, the rotordynamic coefficients are calculated using 
equation (13). The primary measurements are forces, displacements, 
and frequency. The uncertainty in these measurements on the TAMU 
test apparatus are 0.89 N (0.2 lb), 0.0013 mm (0.05 mils), and 0.13 
Hz, respectively. For the four seals tested, the estinoated 
uncertainty in the four rotordynamic coefficients as a percentage of 
the coefficient values was generally less than 5$ for the direct and 
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ROTOR SPCEO <Rp«> 



ROTOR SRCEO <Rp«} 

Fig. 11 Running 8pe«d dependence of direct atlffneas of aeal 1.0 at 
different preaaure ratioa before (above) and after (below) 
noroali ration (low Inlet circumferential velocity In the 
• direction of rotor rotation). See Table H and Fig. 7 for 

deflnltlona of aynbols. 




cross-coupled stiffnesses, from 10 to 20 % for the direct damping, and 
from to several times the value of the cross-coupled damping. 
Since the uncertainties in the cross-coupled damping values were so 
high, and since the cross-coupled-damping forces are of minor 
significance, compared to the other damping and stiffness forces, 
comparisons of the cross-coupled-daraping coefficients have been 
omitted from this report. 

Selection of report data. For each of the four seals tested, there 
were 125 test points for leakage, direct and cross-coupled stiffness, 
and direct damping at the 74.6 Hz shake frequency. Generally, a 
ranking of the three independent variables of the test apparatus in 
order of the relative effect on the normalized rotordynamic 
coefficients of a seal is: inlet circumferential velocity, pressure 
ratio, running speed. Figures in the next two chapters show the 
dependence of leakage and normalized rotordynamic coefficients on 
seal taper for inlet swirl conditions of Table 4. Fig. A1 - Al6 in 
the Appendix show leakage and normalized coefficients at 3000 and 
16000 rpm as a function of inlet circumferential velocity at the five 
applied pressure ratios. Generally, solid lines in a figure 
represent experimental results, and broken lines represent the 
predictions of Nelson's analysis. These figures will be used in a 
comparison of the effect of taper on seal performance, and in an 
evaluation of Nelson's analysis. 


CHAPTER VI 


TEST RESULTS; RELATIVE PERFORMANCE OF SEALS 

This evaluation of the effect on seal performance of varying the 
seal taper angle requires frequent use of the information in Table 4 
and Fig. 7-10. Comparisons of the leakage, direct stiffness, 
cross-coupled stiffness, direct damping, and stability of the four 
seals follow. 

Leakage . The flow rate of air through each seal was measured with 
either a vortex flowmeter or a turbine flowmeter located in the 
piping upstream of the test section (see Fig. i*). The specifications 
on the two meters state ranges of 37 to 53^ and 3 to 60 acfm, 
respectively. Unfortunately, the minimum capability of the vortex 
meter was understated and some erroneous data were recorded, 
predominantly for seal 2.5. Leakage data which were clearly in error 
(e.g. negative leakage) have been eliminated from the figures 
referred to below. 

A comparison of the leakage of the four seals reveals, as 
expected, that for equal exit clearances, greater leakage occurs for 
greater inlet clearances. Data from tests at 3000 and 16000 rpm are 
are shown in Fig. 12 and 13 and in Fig. A1 - A4 in the Appendix. 
Cienerally, the leakage of seals 1.5, 2.0, and 2.5 exceeded that of 
seal 1.0 by from 15 to 25%, 25 to 40$, and 45 to 60$, respectively. 
Direct stiffness. An interesting comparison of the dimensionless 










direct stiffness of the four seals toay be developed from Fig. 14 and 
15 and Fig. A5 - a 6 in the Appendix. The experimental stiffness 
values of seals 1.5 and 2.0 are close throughout the range of inlet 
circumferential velocities. Seal 2.0 appears stiffer at low pressure 
ratios, while seal 1.5 appears stiffer at higher pressure ratios. 
The stiffness of seal 2.5 is clearly less than that of both seals 1.5 
and 2.0 throughout the range of velocities. Seal 2.5 is about 10 to 
40% less stiff than seal 2.0. At a running speed of 3000 rpm, seals 
1.5 and 2.0 are substantially stiffer than seal 1.0 for all 
circumferential velocities. The margin of difference decreases from 
two or three times stiffer to about 20 % stiffer as the 
circumferential velocity increases from negative (opposing rotor 
rotation) to positive. However, at 16000 rpm, seals 1.5 and 2.0 are 
less than 5 % stiffer than seal 1.0 for high pressure, high- 
circumferential-velocity tests. For zero and negative 
circumferential velocities, seals 1.5 and 2.0 are two to three times 
stiffer than seal 1.0 at 16000 rpm. 

Cross-coupled stiffness. Fig. 16, 17, and A9 - A12 in the Appendix 
reveal that seal 1.0 had the greatest and seal 2.5 the lowest cross- 
coupled stiffness magnitudes throughout the range of inlet 
circumferential velocities tested at both 3000 and 16000 rpm running 
speeds. There was little difference in the cross-coupled stiffness 
of seals 1.5 and 2.0 throughout the tests. Furthermore, the cross- 
coupled stiffnesses at zero inlet circumferential velocity were 
positive for seals 1.0 and 2.0 and negative for seals 1.5 and 2.5. 





Pressure: 8.25 bar Pressure: 8.25 bar 

Speed: 3000 rpm Speed: 16000 rpm 
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Fig, IS Experimental normalized direct atlffneases of the teat seals at 8.25 bar pressure. 
See Table 4 for definitions of Inlet circumferential velocity symbols. 


Pressure: 3.08 bar Pressure: 3.08 bar 

Speed: 3000 rpm Speed: 16000 rpm 
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Fig. 16 Experimental normalized cross-coupled stiffnesses of the test seals at 3.08 bar 
pressure. See Table 4 for definitions of Inlet circumferential velocity symbols. 



Pressure; 8.25 bar Pressure: 8.25 bar 

Speed: 3000 rpm Speed: 16000 rpm 
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Fig, 17 Experimental normalized cross-coupled stiffnesses of the test seals at 8,25 bar 

pressure. See Table 4 for definitions of inlet circumferential velocity symbols. 


This data will be discussed further in a stability analysis following 
the comparison of experimental data. 

Direct damping. Fig. 18, 19, and A1 3 “ Al6 in the Appendix show a 
general decrease in the direct damping as the taper (ratio of inlet 
to exit clearance) increases for seals of equal exit clearance. The 
clearest downward shift occurs as the taper is increased from 1 .0 to 
1.5, while seals 1.5 and 2.0 again appear only marginally different. 
Noting that the difference between the damping of seals 1.0 and 1.5 
is significantly less at 16000 rpm, high positive inlet 
circumferential velocity, and high pressures, remember that the 
estimated uncertainty in the experimental damping values was 
generally from 10 to 20 % (see the Section on relative uncertainty). 
Stability analysis. One further parameter of comparison among the 
test seals is the dimensionless whirl frequency ratio. To understand 
the value of this parameter, consider a rotor in a circular orbit of 
amplitude A and frequency 0 (Fig. 20). The X and Y components of 
force in the seal model of equation (3) may be resolved into radial 
and tangential forces 

Fp - Fx cos u)t + Fy sin cot 
Ft “ -Fx sin cot + Fy cos cot . 

Expressing the rotor motion as 

X ■ A cos cot X » “Aco sin (ut 

y •« A sin ut Y - Au cos ut 

and using equation (3), the resultant radial and tangential forces are 
illustrated in the figure and are defined by 


Pressure: 3.08 bar Pressure: 3.00 bar 

Speed: 3000 rpm Speed: 16000 rpm 
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Fig. 18 Experimental normalized direct damping of the teat seals at 3.08 bar pressure. 
See Table 4 for definitions of Inlet circumferential velocity symbols. 


Pressure: 8.25 bar Pressure: 0.25 ba 

Speed: 3000 rpm Speed: 16000 rpm 
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Fig. 19 Experimental normalized direct damping of the teat seals at 8.25 bar oressure. 
See Table 4 for definitions of inlet circumferential velocity symbols. 
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precessing seal rotor. 


If F^/A is a positive quantity, the tangential force is destabilizing 
since it supports the whirling motion of the rotor. Conversely, if 
F^/A. is negative, it opposes the whirling motion of the rotor, and is 
therefore stabilizing. The whirl frequency ratio is defined by 
Whirl frequency ratio - k / CQ . 

From the above discussion, if the whirl ratio is less than one, the 
tangential force on the rotor is stabilizing. A minimum value of the 
whirl frequency ratio is optimum for stability. 

Fig. 21 - 24 show the whirl-frequency ratios at a running speed 
of 16000 rpm and a shake frequency of 74.6 Hz of the seals repiorted 
on here. From the lowest whirl frequency ratios to the highest, the 
seals rank as follows: seal 2.5, seal 1.5, seal 2.0, and seal 1.0. 



INLET circumferential VELOC 1 TY (•/•«€> 

Fig. 21 Whirl frequency ratios of seal 1.0. Five pressure ratios, 
rotor speed - 16000 rpm, shake frequency - 7**. 6 Hz 


a 

2.7 
2 . 4 


1 



INLET CIRCUMFERENTIAL VELOCITY 

Fig. 22 Whirl frequency ratios of seal 1.5. Five pressure ratios 
rotor speed - 16000 rpm, shake frequency - 7^1 *6 Hz 
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inlet CIRCUMPERENTIAL VELOCITY (m/sac) 

23 Whirl frequency ratios of seal 2.0. Five pressure ratios 
rotor speed - 16000 rpa. shake frequency • 7**.6 Hz 


inlet circumferential velocity Wsac} 


24 Whirl frequency ratios of seal 2.5. Five pressure ratios 
rotor speed - 16000 rpsL. shake frequency • 74.6 Hz 




CHAPTER VII 


TEST RESULTS; COMPARISON TO THEORETICAL PREDICTIONS 

In this chapter, the experimental results from the tests of four 
smooth-rotor, smooth-stator, tapered seals are compared to the 
predictions of an analysis developed by Nelson [10]. The seals 
tested are described in Tables 1-3* Table and Fig. 7 ~ 10 define 
the symbols used in the figures. 

Dependence on inlet loss coefficients. As noted previously. Nelson’s 
analysis has underpredicted the direct, stiffness of tested seals in 
previous comparisons. Nelson et al. [11] have suggested that 
inadequate modeling of the seal inlet losses could account for the 
underpredictions. Of the rotordynamic coefficients, direct stiffness 
is the most sensitive to changes in the entrance loss coefficient, k, 
of equations (8); and (9). This statement is supported by the results 
of Fig. 25. In this figure, the broken lines represent the 
normalized rotordynamic coefficient values determined during a 

particular test of seal 1.0. The solid lines represent the values 

predicted by Nelson's Moody-based analysis [10] as the entrance loss 
varies from 0.1 to 0.9. Recall that in the entrance loss 
relationship of equation (8), 

_ 2 Y/d-Y) 

Po(0) - {1 [(Y-1 )(k+1)Mo (0)] / 2} (8) 

PqCO) is the seal entrance/reservoir pressure ratio and Mo(0) is the 
entrance Mach number. In the past, these values, and therefore the 



DAMPING (imac)- STIFFNESS (noodJmi 


entrance loss, k 


ENTRANCE LOSS. 


Dependence of stiffness (above) and damping (below) on 
entrance loss coefficient. Curves 1 : direct stiffness 
(Kxx) and damping (Cxx). Curves 2 : cross coupled (Kxy 
and Cxy). Broken lines represent experimental data from a 
test of seal 1.0: rotor speed - 16000 rpm, inlet 
circumferential velocity - 45 m/sec, pressure = 8.25 bar. 





value of k, have been determined from pressure measurements at the 
first pressure tap in the inlet end of each seal. In Fig. 26, a 
pressure gradient plot is presented for the tests of seal 1.0 which 
is typical of all seals tested. The shape of the curves show that 
some of the entrance loss occurs after the initial sharp drop in 
pressure at the seal inlet. The first two data points on the curves 
represent the reservoir pressure for the tests. The remaining points 
represent the pressures measured at taps along the seal. Nelson's 
analysis assumes that the entrance loss occurs at the entrance. In 
ein attempt to represent the entrance losses which actually occur in 
the first centimeter of the seal, a curve fit of the form 

e 

p(x) - P(a - x) (16) 
was applied • on the data points which follow the inlet loss. In 
equation (16), x is the axial location of the absolute pressure p(x), 
and P* a, and 6 are obtained by a least squares curve fit. The value 
of p(0) is then used as the seal entrance pressure in Po(0) of 
equation (8) to determine the entrance loss k. Fig. 27 shows an 
example of the entrance losses thus obtained. The remaining entrance 
loss figures may be found in the Appendix, Fig. A17 - A26. Before 
leaving Fig. 26, it should be noted that five curves are plotted, one 
for each test rotor speed. This figure illustrates the typically 
negligible effect of rotor speed on the pressure distribution in the 
seal. Subsequent pressure plots in this report will include data 
from one test speed only. 
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Fig. 27 Entrance los<a for seal 1.0, high velocity swirl in the 
direction of rotor rotation. 
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Fig. 28 Direct stiffness, comparison of Nelson's analyses (low 
velocity swirl in the direction of rotor rotation). 

Solid - experimental; dash - Moody; dash-dot-dash ^ Hlrs 


Comparison of Nelson’s analyses. In the comparison of his Hirs-based 
and Moody*-based analyses, Nelson [10] used a water seal for input 
parameters. He found that the use of Moody's equation resulted in 
higher direct and cross-coupled stiffness predictions for "rough" 
seals (e/2HQ ■ 0.05), but that the differences were small for 

relative roughness near zero. There was little difference in the 
direct damping predictions of the two analyses. 

Fig. 28 - 30 show a comparison of the predictions of Nelson's 
analyses using input data from the tests of seal 1.0. In these 
figures, solid lines represent experimental data, dash lines are the 
predictions of Nelson's Moody-based analysis, and dash-dot-dash lines 
are the predictions of Nelson's Hirs-based analysis. Though the 
measured relative roughnesses of the rotor and stator of seal 1.0 
were only 0.0003^ and 0.00050, respectively, the direct stiffness 
predicted by Nelson's Moody-based analysis was from 10 to 25!t higher 
than predicted by the Hirs-based analysis. The cross-coupled 
stiffness and direct damping predictions differed by less than 10?. 
Static results comparison. Fig. 31 - 38 and Fig. A1 - A16 in the 

Appendix are useful in the evaluation of Nelson's analysis which uses 
Moody's equation. In the comparison plots, the theoretical 
predictions and experimental results are represented by broken and 
solid lines, respectively. The numbers on the plotted lines refer to 
the inlet conditions listed in Table Comparisons of leakage 

predictions to experimental results may be found in Fig. 3''. 32, and 
A1 - A4. Predictions of increasing leakage with increasing inlet 
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Speed: 3000 rpm Speed: 16000 rom 



Fig. 32 Experimental and theoretical mass flow rates of air through the test seals at 8.25 bar 
pressure. See Table 4 for definitions of circumferential velocity symbols. Solid 
lines - experiment; broken lines - theory. 




clearance agree with experimental result?. The prediction? were 
within 5 % of the measured mass flow rate in some cases, and always 
within 

Plots of the theoretical and experimental pressure gradient in 
each seal are in Fig. A27 - A37 in the Appendix. The plots show good 
agreement between theory and experiment. Most of the predicted 
pressure gradients for the seals have larger pressure drops at the 
inlet than do the experimental gradients because the entrance loss in 
Nelson's model takes place at the inlet, while experimentally it 
appears to take place in the first centimeter of the seal. Also, the 
slopes of the experimental gradients are steeper than the predicted 
gradient slopes. One explanation is that the experimental gradients 
were measured along the stator, which was rougher than the rotor in 
all cases, while the theoretical gradients are developed from the 
combined effect of the rotor and stator roughnesses. A final 
observation may be made concerning the erratic shape of the 
experimental pressure gradients of seals 1.5 and 2.0. While the 
sixteen pressure taps of seals 1.0 and 2.5 were aligned, as in Fig. 
6, the arrangement was different for seals 1.5 and 2.0. In these 
seals, which were fabricated first, the first eight pressure taps 
were aligned as in Fig. 6, and the last eight were aligned 90® 
counter-clockwise from the first eight. The sharp drop in pressure 
after the eighth pressure tap in the experimental pressure gradients 
of seal 1.5 and the similar shape for seal 2.0 are thought to be a 
result of this arrangement. The selection of data used in the 


pressure curve fits described in the section on entrance losses was 
made more difficult by these phenomena. 

Dynamic results comparison. The dynamic results compared include the 
normalized experimental and theoretical direct and cross-coupled 
stiffness and direct damping coefficients. A cross-coupled damping 
comparison has been omitted because of the uncertainty present in the 
experimental values (see the Relative uncertainty section). Of the 
remaining three coefficients, the direct stiffness comparison will be 
presented first, and the cross-coupled stiffness comparison will be 
last. 

Direct stiffness. Fig. 33, 3^. and A5 - A8 in the Appendix show the 
theoretical and experimental direct stiffnesses of the four seals 
tested. Fig. 33 and 3^ show the theoretical and experimental 
dependence of direct stiffness on seal taper. The data plotted are 
for inlet swirl configurations 3 and 5 of Table Nelson's analysis 
predicts significantly higher direct stiffnesses for seals 1 .5, 2.0, 
and 2.5 than it predicts for seal 1.0. Furthermore, the analysis 
predicts that seal 1 .5 is the stiffest seal and does not predict that 
seal 2.0 is stiffer than seal 2.5 for high-speed inlet swirl in the 
direction of rotor rotation. Experimentally, little difference was 
present in the stiffness of seals 1.0 and 1.5 at 16000 rpm, seal 2.0 
was generally the stiffest seal, and seals 1.5 and 2.0 were always 
stiffer than seal 2.5. 

The greatest difference between theoretical and experimental 
direct stiffness coefficients for a seal is seen in Fig. A5, where 
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agreement for seal 1 .0 Is best for inlet air swirl opposing the rotor 
rotation direction. The experimental direct stiffness was more 
speed- and swirl-velocity dependent than predicted. For swirl in the 
direction of rotor rotation, theory underpredicts the experimental 
values by 30 to 60$ at the 16000 rpm rotor speed. The direct 
stiffness predictions for the remaining seals were generally within 
30$ of the experimental, values, with a number of predictions within 
10$. The greatest disagreement is in the non-prerotated case for 
each seal. Experimentally, the direct stiffness for a given inlet 
pressure was lowest for this case, while theory predicts a higher 
value for the non-prerotated case than for the two low speed swirl 
cases. 

Direct damping. In the comparison of the normalized theoretical and 
experimental direct damping data (Fig. 35, 36, and A13 - Al6 in the 
Appendix), recall that the estimated uncertainty was generally from 
10 to 20$ of the experimental direct damping. The predictions of 
Nelson's analysis were generally within 30$ of the experimentally 
determined direct damping of seals 1.5, 2.0, and 2.5. The same is 
true of seal 1.0 at the three higher inlet pressures of Table 
However, at the two lower pressures, theory underpredicts experiment 
by up to 60$. A general trend of decreasing damping with increasing 
taper is present in both experimental and theoretical results. 
Cross-coupled stiffness. Fig. 37 and 38 show that the experimental 
cross-coupled stiffness of seal 1 .5 was 30 to 50$ less than that of 
seal 1.0 for inlet swirl configuration 5 of Table 4. Nelson's 


Pressure: 3.08 bar Pressure: 3.00 bar 

Speed: 3000 rpm Speed: 16000 rpm 
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Fig. 35 Experimental and theoretical normalized direct damping of the teat seals at 3.08 bar 
pressure. See Table 4 for definitions of circumferential velocity symbols. Solid 
lines - exfperlment; broken lines - theory. 
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Fig. 36 Experimental and theoretical normalized direct damping of the test seals at 8.25 bar 
pressure. See Table 4 for definitions of circumferential velocity symbols. Solid 
lines - experiment; broken lines - theory. 
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analysis predicts little difference between the cross-coupled 
stiffnesses of seals 1.0 and 1.5. In addition, the experimental 
cross-coupled stiffness values of seal 2.0 were higher than those of 
seals 1.5 and 2.5, but Nelson’s analysis predicts that seal 2.0 has 
lower values than seals 1.5 and 2.5. 

For a given seal. Nelson's analysis predicts an almost linear 
relationship between cross-coupled stiffness and the circumferential 
velocity of the inlet air (Fig. A9 - A12 in the Appendix). 
Experimentally, the relationship is near linear through the low speed 
swirl in the direction of - rotor rotation. The linearity is not 
present when the inlet swirl speed is further increased. As a 
result, the normalized cross-coupled stiffness predictions are as 
much as seven times (seal 2.5) greater than the experimental values 
at the highest circumferential velocity of the inlet air. Agreement 
is best when the magnitude of the inlet swirl velocity is less than 
eighty meters per second. 

Stability analysis. The whirl frequency ratio, a measure of seal 
stability, was defined in an earlier section. Fig's. 39 - ^2 present 
the whirl ratios predicted by Nelson's analysis for the seals at the 
16000 rpm rotor speed. Five curves are present in each figure, one 
for each inlet pressure applied (see Table ^) . Fig's. 21 - 24 
provide the corresponding plots of experimental whirl ratios. For 
inlet swirl in the direction of rotor rotation (positive 
circumferential velocity of the inlet air), the most frequent 
situation in practice. Nelson predicts less stability for each seal 


WHIRL FREQUENCY RArTlO ^ WHIRL FREQUENCY RATIO 













than was found experimentally. This is a result of cross-coupled 
stiffness overpredictions and direct damping underpredictions. The 
difference in predicted whirl frequency ratios among the seals are so 
small at a given inlet circumferential velocity that stability 
comparisons are inconclusive. 



CHAPTER VIII 


CONCLUSIONS 

An air seal teat apparatus at Texas A & M University has been 
modified to enable tests at rotor speeds up to I6OOO rpm. A 
previously established method of determining seal rotordynamic 
coefficients has been used to obtain consistent, repeatable stiffness 
and direct damping results. 

The experimental and theoretical results of the preceding 
section support the following conclusions: 

(a) Four smooth-rotor/smooth-stator seals with equal exit 
clearances were tested. The inlet to exit clearance ratios (tapers) 
of the seals were 1.0, 1.5, 2.0, and 2.5. The leakage through seals 
1.5, 2.0, and 2.5 were greater than through seal 1.0 by 15 to 25?, 25 
to M0$, and i»5 to 60?, respectively. 

(b) Throughout the tests, the differences in the rotordynamic 
behavior of seals 1.5 and 2.0 were small. 

(c) Generally, seals 1.5 and 2.0 had clearly higher normalized 
experimental direct stiffness coefficients than seals 1.0 and 2.5. 
At a rotor speed of 16000 rpm, seal 1.0 was sometimes stiffer than 
seal 2.5. 

(d) The direct stiffness of seal 1.0 is speed dependent for 
inlet air swirl in the direction of rotor rotation. At rotor speeds 
in excess of I6OOO rpm, seal 1.0 might prove stiffer than seals 1.5 
and 2.0 for inlet swirl in the direction of rotor rotation. 


(e) A comparison of the experimental cross-coupled stiffness 
and direct damping of the seals through the whirl frequency ratio (a 
nondimensional stability indicator) leads to the following seal 
ranking, in order of decreasing stability: seal 2.5, seal 1.5, seal 
2.0, and seal 1.0. 

(f) For rotor speeds up to 16000 rpm, the effects of speed 
Increases on the rotordynaroic coefficients of the smooth- 
rotor/smooth-stator seals tested are mainly due to the seal clearance 
reduction as the rotor grows. 

(g) The evaluation of seal entrance losses is critical when 
using Nelson's analysis to predict seal performance. Aside from 
direct test data, no readily available source of information in this 
area is known. 

(h) Although the experimental direct stiffness of seal 1.0 was 
always lower than that of seals 1.5 and 2.0, the uncommon speed 
dependence of seal 1 .0 made the difference modest for inlet air swirl 
in the direction of rotor rotation at a rotor speed of 16000 rpm. 
Nelson's analysis indicates that seals 1.5 and 2.0 are significantly 
stiffen than seal 1 .0 at all speeds and inlet swirls tested. In 
addition. Nelson's analysis predicts the normalized direct stiffness 
of seal 2.5 to be slightly higher than that of seal 2.0 for high 
speed inlet air swirl in the direction of rotor rotation, but seal 
2.0 was always clearly stiffen than seal 2.5 experimentally. 

(i) The normalized direct damping predictions of Nelson's 
analysis are generally within 30% of the experimental results. A 


trend of decreasing damping with increasing taper is present in both 
experimental and theoretical results. 

(J) Nelson's analysis predicts little difference between the 
cross-coupled stiffnesses of seals 1.0 and 1.5. Experimentally, the 
normalized cross-coupled stiffness of seal 1.5 was 30 to 50% less 
than that of seal 1 .0. Furthermore, Nelson's analysis predicts a 
linear relationship between cross-coupled stiffness euid the 
circumferential velocity of the inlet air, which is not present 
experimentally. Cross-coupled stiffness predictions are as much as 
seven times greater than experimental values for high inlet 


circumferential air velocities. 
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inlet circumferential velocity (m/sac) 



inlet circumferential velocity («/eec) 


Fig. A2 Maas flow rate of air through seal 1.5 at 3000 rpm (above) 
and 16000 rpm (below). See Table ^4 for pressure symbol 
definitions. 
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inlet circumferential velocity (m/soc) 
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Fig. A14 Normalized direct damping of seal 1.5 at 3000 rpm (above) 
and 16000 rpm (below). See Table 4 for pressure symbol 
definitions. 
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Fig« Al6 Morasllzed direct demping of seal 2.5 at 3000 rpo (above) 
and 16000 rpo (below). See Table 4 for pressure symbol 
definitions. 
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Fig. A20 Entr^ce losses of seal 1.5, inlet air swirls 3 (top) and 

(botton) of Table 4. Curves are numbered by rotor speeds 
of Table 4. 
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Fig. A26 Entrance losses of seal 2.5, inlet air swirl 5. Curves 
numbered by rotor speeds of Table 4 . 
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Fig. A27 Pressure gradients of seal 1.0, inlet swirl 1. Solid 
lines - experimental data; broken lines - theory. 
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Fig. A36 Pressure gradients of seal 2.5, Inlet air swirls 3 (top) and 
4 (bottom). Solid lines - experiment; broken lines - theory 
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Fig. A37 Pressure gradients of seal 2.5, inlet air swirl 5. 

Solid lines - experiment; broken lines - theory. 
See Table ^ for definitions of symbols. 
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